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ABSTRACT 

After reviewing the thermodynamics and critical phenomena associated with 
AdS black holes carrying multiple R-charges in various dimensions, we do a 
Bragg- Williams like analysis of the systems around its critical points. This 
leads us to propose an effective potential governing the equilibrium properties 
of the boundary gauge theory. We also study certain non-equilibrium phe- 
nomena associated with these gauge theories. In particular, we compute the 
conductivities and diffusion coefficients for theories with multiple R-charges 
in four, three and six dimensions. 
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1 Introduction 



Recently there appear stimulating studies concerning spinning near-extremal D3 branes 
[1-9]. There are at least two reasons that motivate these activities. The more important 
one is, perhaps, the fact that this provides a simple example in the framework of gauge- 
gravity duality where different aspects can be studied in presence of non-zero charge as 
well as chemical potential. Indeed, the ten dimensional gravitational background corre- 
sponding to near extremal spinning D3 brane provides gravity dual of a four dimensional 
Af = 4 finite temparature super Yang-Mills theory. The number of independent angular 
momenta of the spinning D3 brane is the rank of the isometry group 5*0(6) of the space 
tranverse to the brane, which upon dimensional reduction corresponds to charged black 
hole with three U(l) charges. In the gauge theory, this isometry group represents the 
R-symmetry group. Three independent chemical potentials can therfore be introduced 
which couple to the three independent U(l) charges. By exploiting the gauge-gravity 
correspondence, several recent works explored the behaviour of such gauge theories at 
strong coupling in its simplest setting. These are the cases where only one of the chemical 
potentials is non-zero. 

Another reason that makes the study of rotating D3 brane interesting is the fact 
that spinning D3 brane has an upper critical limit of angular momentum beyond which 
thermodynamic instability sets in. Though the final configuration to which this instability 
leads is yet to be understood, it is conjectured that when angular momentum reaches 
this critical value, the system undergoes a phase transition. A careful inspection of the 
bahaviour of various thermodynamic quantities suggests that it is a second order phase 
transition. In particular, near this critical point, divergence of correlation length shows up 
in the behaviour of thermodynamic quantities such as specific heat and non-equilibrium 
quantities such as, transport coefficients. One expects that, in the vicinity of the second 
order phase transition, certain physical properties of the system can be extracted without 
the knowledge of their microscopic behaviours. 

With these motivations, in this paper, we study various aspects of the R-charged 
black holes near its critical point. So far these studies were carried out for the simplest 
variant of this model, namely, D3 branes with a single R-charge. The purpose of this article 
is to extend these stuides to its full generality. This leads to more intricate phase structures 
of the systems along with interesting surprises. So far the equlibrium thermodynamics is 
concerned, we carry out a Bragg- Willinas like analysis in the bulk. More specifically, by 
identifying the horizon radius as the order parameter, we write down the Bragg- Williams 
[10] potential for the R-charged black hole. We also extend our analysis for charged black 
holes in four and seven dimensions. They appear as the near horizon limits of non-extremal 
spinning M2 and M5 branes. Subsequently, by exploiting the gauge-gravity duality, we 
illustrate with an explicit example, how this leads to a proposal of effective potential 
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for the gauge theory. As far as the non-equilibrium phenomena are concerned, we study 
the transport coefficients of the gauge thery duals using gauge-gravity correspondence. 
Besides giving a general precription to solve perturbed gravity equations for multiple R- 
charged black holes in four, five and seven dimensions, we explicitly compute the diffusion 
coefficients for all these cases using the Green-Kubo formula. 

The paper is structured as follows. In the next section, we reconsider various R- 
charged black holes, in the grand canonical enseble, with a focus on their thermodynamic 
properties [2]. On one hand, this section allows us to set our notations and conventions and 
on the other hand, it helps us to scan the complicated phase structure associated with 
the black holes with multiple R-charges in various dimensions. Next, in section 3, we 
consider the black holes near their critical temperatures. By identifying horizon radius as 
the order parameter, we construct the Bragg- Williams like potentials for the black holes. 
This further leads us to propose an effective potential for the boundary gauge theory 
with charge density as the order parameter. We illustarate this explicitly for the four 
dimensional gauge theory with a single non-zero chemical potential. Section 4 is devoted 
to a study of certain non-equilibrium phenomena associated with these bulk geometries 
and its imprint on the boundary theories. In particular, we compute transport coefficients 
using Green-Kubo formula. It turns out that for multiple R-charges, the calculations are 
quite involved. We, therofore, in the main text restrict ourselves to the evaluation of 
the R-charge conductivity for black holes when two charges are non-zero deferring rest 
of the cases to the appendix. The appendix also includes an algorithm to compute the 
conductivities in various dimensions and for multiple charges. Finally, we conclude the 
paper in section 5 with a discussion on our results. 

2 Black Holes in STU model 

In this section we review thermodynamics of the R-charged black hole in five, four and 
seven dimensions. We start with the black holes in five dimensions. These are the solutions 
of equations of motion of D = 5, M = 2 gauged supergravity obtained by compactification 
of ten dimensional IIB supergravity on S 5 . These are known as STU black holes and were 
found in [11-13]. The effective action is given by 

(2.1) 

where / represents the scale associated with the cosmological constant. In addition to 
the metric, we have three scalar fields X\ i = 1,2,3 which are constrained through the 
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relation X 1 X 2 X 3 = 1. V is a potential involving the scalar fields given by 

V = 2 t,X-r ( 2 - 2 ) 
1 

Gij represents metric of the moduli space parametrized by the scalars and is given by 

Gv = ±dbe[(x 1 )- 2 ,(x 1 )-*,(x 1 )- 2 ]. (2.3) 

F l are the field strengths associated with three abelian gauge fields A 1 corresponding to 
three U(l) in the cartan of the R-symmetry group. 

As shown in [11], this effective action (2.1) admits asymptotically AdS black hole 
solutions with three U(l) charges. These solutions can be written down as 

2/2 2 

ds 2 = -H-^fdt 2 + H^dr 2 + H^(dx 2 + dy 2 + dz 2 ), (2.4) 
V r l l z 

where 

/ = <! + > + > + £)-£ (2-5) 
and the harmonic functions are given by 



H = H 1 H 2 H 3 =(1 + ^)(1 + ^)(1 + ^). (2.6) 

In the above solutions, the parameter tq is related to the mass or energy density of the 
black hole while qi,q2, <?3 are related to the charges of the black hole. 



2.1 Black hole thermodynamics and instabilities in five dimen- 
sions 

In this subsection we briefly describe the thermodynamics and the critical phenomena 
associated with the five dimensional R-charged black hole given by (2.4,2.5,2.6). From 
now, instead of using r as a parameter, we will use the radius of the outer horizon r + 
to parametrize the solution. r + can be obtained from the largest root of the equation 
f(r+) = 0, where f(r) is given in (2.5). Furthermore, we will scale all the dimensionful 
parameters with AdS length scale / e.g. we introduce the dimensionless parameter f = 
r + / 1 and T — IT for horizon radius and Hawking temperature respectively. Similarly the 
dimensionless charge parameters will be denoted by qi = qijl 2 . In terms of these scaled 
parameters, horizon area becomes proportional to 

A = v /(r2 + g 1 )(f2 + g 2 )(f2 + g 3 ). (2.7) 
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In addition, we will set Newton's constant G — 1. G and I will be resurrected whenever 
required. In what follows we will consider the system in the grand canonical ensemble 
with temparature T and chemical potentials /ij, i — 1, 2, 3 as external parameters. 

The scaled Hawking temperature turns out to be 

2f 6 + f 4 (gi + q 2 + q 3 ) - q x q 2 q 3 
T ~ 2^4 • (2 ' 8) 

Other thermodynamic quantities such as entropy, energy and pressure are expressed as 

where the above quantities are rendered dimensionless by multiplying appropriate power 
of I. The components of the charge density and chemical potential are given by 

- ,,, 1. (2.10) 



IQnr r(r 2 + q^) 

Since we are analysing black holes with flat horizon, the horizon has infinite volume. 
Consequently, only the themodynamic densities are finite. What we have written above 
actually represent those thermodynamic densities. 

Unless charges satisfy certain constraints, these black holes undergo a local insta- 
bility [2,3,14]. While at high temperature, black holes remain stable, once we reduce the 
temperature down to a critical value, the specific heat and suceptibility diverge. In order 
to see this, let us compute those quantities. The specific heat associated with the black 
holes has the following form 



C = (f-p) = (2f 6 + f 4 (g x + q 2 + q 3 ) - q.Ms) X 

3f 6 - f 4 (gi + q 2 + q 3 ) - f 2 (q x q 2 + q 2 q 3 + Mi) + MiMs 
4A(2f 6 - f 4 (g! + q 2 + q 3 ) + q^qs) 

The susceptibility, for the present case is a 3 x 3 symmetric matrix, defined as 



(2.11) 



xa= Gf 1 )- , with k,l^j. (2.12) 

Two of its components are given by 

_ 2f 8 + 5f 6 gi - f 4 qj - f 6 q 2 - f A q x q 2 - f 6 q 3 - f 4 <M 3 - 3f 2 gig 2 ?3 + QiMs 
XU 167r(2f 6 - f 4 (g x + q 2 + q 3 ) + qiMz) ' (0 , ^ 

r^=- -4- ,-4- -4- K^- 16 ) 

^[q^q 2 r q 1 + r q 2 - r q 3 - q x q 2 q 3 
Xl2 8vr 2f 6 - f 4 (qi + q 2 + q 3 ) + M2?3 ' 
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The other components can be obtained by cyclically permuting the indices (123). We 
now note that the specific heat and the suceptibility components diverge over the critical 
hypersurface 

2f 6 - f\q 1 + q 2 + q 3 ) + q x q 2 q z = 0. (2.14) 

The thermodynamic quantities for a subclass of black holes where one or more 
charges are turned off can be easily obtained from the above by appropriate substitutions. 
Apart from some isolated cases, which we will discuss in the sequel, all these black holes 
show instabilities. The only difference is that hypersurfaces at which divergences occur 
are different from (2.14) and depend on the nature of the black hole. This is listed below 
for black holes with single, double or three charges : 

2f 2 - qx = 0, for q 1 ^ 0, q 2 = <? 3 = 0, 
2f 2 - q 1 - q 2 = 0, for q~i,q 2 ^ 0, q 3 = 0, 

2 f 4 - q,f 2 - M2 = 0, for q 1 ^0,q 2 = q 3j t 0. (2.15) 

We also note that thermodynamic quantities are everywhere well behaved in the cases 
when charges satisfy certain restrictions. These include q~\ = q 2 , q 3 = and q± = q 2 = q 3 . 
The second one is the standard Reissner-Nordstrom black hole in AdS space. Thero- 
modynamic properties of these black holes were discussed in detail in [3]. The critical 
surfaces given in (2.15) can be expressed as relations between temperature and chemical 
potentials. For example, for one and two charge black holes, they are given respectively 
by 

rp jh 
J- 1 
7T 

j, = {fii + /ff + (fti - £2)3 (fti + + (fii - £2)^1 + ^2)^ 

(2.16) 

For a generic three charged black hole, we have not been able to find the critical surface 
in terms of temperature and chemical potentials in a compact form. For a given set of fii, 
(2.16) gives us the critical temperature at which black hole instability arises. As we do 
not yet know the stable configuration below this temperature, in this paper we consider 
the black holes slightly above the critical temperature. 

As one approaches the critical surface the correlation length diverges. This shows 
up in the divergences of thermodynamic quantities. Near the critical temperature, these 
black holes have universal features and their behaviour do not depend on details of the 
theory. These features are encoded in the critical exponents associated with behaviour of 
various thermodynamic quantities as we approach the critical surface. 
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We begin with the following four static critical exponents. Approaching the crit- 
ical surface with fixed /ij, one finds three critical exponents a,/3, 7 associated with the 
bahaviour of specific heat, charge densities and susceptibilities respectively. 

C ~ (T - f c )- a , p - p c ~ (f - T^, x ~ (f - f c H. (2.17) 

Since there are three C/(l)'s, there are three charge densities and so the number of respec- 
tive exponents $ is also three. However, it turns out that they are all equal and one gets 
only a single (5. Similarly, the susceptibility matrix Xij being symmetric and real, can be 
diagonalized as Xij — Xi$ij an d generically one would expect three 7's one for each of the 
components. But since all the three components are equal we get one 7 only. On the 
other hand, exponent 5 is obtained by approaching the critical line with a trajectory on 
which the temperature is constant. This is defined as 

p-pc (fi-p c )^. (2.18) 

We find for all the black holes which show instability have the same set of critical exponents 

(a,/?, 7 ,*) = (^^2). (2.19) 

These critical exponents are not independent but are related through scaling relations. 
These critical exponents satisfy 

a + 2/3 + 7 = 2, 1 = (3(5-1). (2.20) 



There are two other critical exponents v and rj associated with the behaviour of 
correlation length and correlation function near critical surface. If G(r) is the correlator 
its behaviour near the critical surface is 



G(r)~e-^ at T ^ T c , £ 
G(r) ~ r - d+2 -\ at T = T c . 



T-T c 

T c ) (2.21) 



Assuming additional scaling relations 7 = v(2 — 77) and 2 — a = vd, to be valid we can 
determine v and rj to be equal to - and 1 respectively. However, it is not clear a priory 
whether the scaling relations are valid. One can compute v and i] from correlation of 
scalar modes in the gravity theory. However, if we assume the scaling relations are valid 
then this will belong to class of B-model in classification given in [15]. 

As we will see in the next section, it is possible to construct a Bragg- William en- 
ergy function which reproduces all these critical exponents near the second order phase 
transition surfaces. This function, in turn, will allow us to propose an effective potential 
for the boundary gauge theory via AdS/CFT correspondence. 
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2.2 Black Hole thermodynamics and instabilities in four and 
seven dimensions 



In this subsection, we briefly discuss the thermodynamics and instabilities of R-charged 
black holes in four and seven dimensions. As mentioned earlier, these configurations are 
obtained from spinning Ml and Mb branes of eleven dimension supergravity. First, let 
us focus on the black holes in four dimensions. They can carry atmost four independent 
charges arising form Cartans of SO (8). The black hole solutions are discussed in [12]. For 
our purpose, we require the thermodynamic quantities only. Our notation remains the 
same as in the last subsection except now there will be four charges q iy % = 1,2, 3, 4. For 
convenience, once again we introduce the scaled variables f and cji and the horizon area 
is proportional to 

A = sj{f + q 2 )(f + q 2 )(f + q 2 )(f + q 2 ). (2.22) 
The temperature is given by 

f = [-I+V—LJa (2.23) 

The entropy, energy and pressure are given by 

S = -A , E = —A 2 , P = —^—A 2 . (2.24) 
. The components of charge and chemical potentials are 



J_ /2f / 2q 2 

32vrV f ^' ^-Vf(f + q 2 y 



p* = ^z\l-r A > fii = -2^ A - ( 2 - 25 ) 



The specific heat is given by 

e =( f I) =JKa< < 2 - 26 ' 

where ci, C2, A are given by 

o-4 / -2\-3 , /-2-2-2 , -2-2-2 , -2-2-2 , -2-2-2\- o TT -2 

c 1 = 2r - ( 2^ Qi )r + (q^Qs + + %^<li + IaIi^F ~ 2 H Qi , 

i=l,4 i=l,4 

c 2 = 3f 4 + 2(J2 tV + (tt + tt + tt + tt + lit + tt)r 2 - II t, (2.27) 

i=l,4 i=l,4 
A = 3f 4 - 2( ^ 3 + fe 2-2 + -2-2 + -2-2 + -2-2 + -2-2 + -2 -2 )f 2 _ "Q -2 

i=l,4 i=l,4 

The suceptibility in this case is a four dimentional metrix defined as 

Ki = ilT-)- ' with ( 2 - 28 ) 

VC/ij / T,iA, k ,m,n m 
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Two of the components are 

• dpi \ f5i 



\dfii) T,fi2,p3,p4 327rA' 



Xl2 =lo— J- = ~ oo A ' ^ 2 ' 29 ) 



where 



ft = 3f 4 + 3(2ft 2 - ft 2 - ql - ft 2 )f 3 + (ft 2 ft 2 + gfft 2 + (&q\ - 3(ft 2 ft 2 - ft 2 ft 2 - q 2 qj)) + 3 J[ ft 2 , 
and 



r -4 , o/-2 i -2 -2 -2\-3 , /o-2-2 -2-2 -2-2 -2-2 -2-2 -2-2\-2 , TT - 

5 2 = r + 2(ft + ft - ft - ftF + (3ftft - ft ft - QiQs ~ - ft ft ~ 9 2 ft) r + [[ 5. 



■2 

j=l,4 



Other components of suceptibility can be obtained by cyclically permuting the indices 
(1234). We note that specific heat and suceptibility diverge at A = 0. 

Apart from some isolated cases (discussed below), all black holes suffer from insta- 
bilities. When some of the charges are either equal or zero, the critical surfaces change 
as can be seen directly by sustituating charges in the general expressions above. We list 
down various such critical surfaces below when some of the charges are identically zero: 

3f - 2ft 2 = 0, for ft ^ 0, ft = ft = ft = 0, 

3f 2 - 2f (ft 2 + q 2 2 ) + q\q\ = 0, for ft, ft ^ 0, ft = ft = 0. (2.30) 

Furthermore, unlike the black holes in five dimensions, these black holes show instability 
even when two of their charges are equal. In particuler, ft = ft 7^ and with ft = ft = 0, 
the singularities appear at 

3f - 2ft 2 = 0. (2.31) 
For black holes with three charges, the critical surfaces are 

3f 2 - 2f (ft 2 + ft 2 + ft 2 ) + ft 2 ft 2 + ft 2 ft 2 + q%q% = for ft, ft, ft ^ 0, ft = 0, 

3f - qf - 2ft 2 = for ft = ft ^ 0, ft ^ 0, ft = 0. (2.32) 

We further note that, for three non-zero charges, black holes are thermodynamically stable 
only when ft = ft = ft. Though for generic values of four charges, instabilty appears 
at A = 0, when some of the charges are equal location of thermodynamic singularities 
change. This is given in the following list. 

3f 3 - f 2 (2ft 2 + 2ft 2 + ft 2 ) + fq 2 x ql + ft^ft 2 = for ft + ft, ft = ft, 
3f 3 - f (ft 2 + ft 2 ) - q\q\ = for ft = ft, ft = ft, 

3f 2 - 2fft 2 - q\q\ = for ft = ft = ft, ft + 0. (2.33) 



8 



However, when all the charges are equal, thermodynamic instabilities are absent. 

Approaching the critical surfaces in different ways, it is possible to find all the critical 
exponents as in the previous subsection. In all the cases it turns out that the exponents 
are same as in (2.19). 

In the rest of this subsection, we consider the thermodynamic properties of black 
holes in seven dimensions. Since those black holes correspond to S 4 reduction of spin- 
ning M5 branes in eleven dimensional supergravity, the Cartan consists of U(l) x U(l). 
Therefore, a general R-charged black hole can carry only two independent charges. The 
general form of the solutions are given in [12]. We only list down the thermodynamic 
variables. The horizon area is proportional to 

A = yV + g^ + gl), (2.34) 

The temperature is 

_ = 3f 8 + f\ql + q 2 2 ) - q\q\ 

27Tf 3 A ' 1 ' ' 

Entropy, energy and pressure are 

The components of the charge and chemical potential are 

Pi = , ft = A A (2.37) 

8nr r(r 4 + qf) 

The specific heat is 

- f(5f 8 - f 4 q\ - f 4 q 2 2 - 3g 1 2 g|)(3f 8 + f 4 qj + f A q% - q\q\) 

3f 8 _ f ±q\ - f 4 q\ - q\q\ ' 1 ' 

The suceptibility is a two dimensional symmetric matrix with entries 

f(3f 12 + 13f 8 g 2 — f 8 g| — 4f 4 qf + f 4 q\ql — ^qfq^) 



Xn = 



37r(3r 8 — r 4 q\ — r 4 q\ — q\ql) 



qiq 2 r{5qiq2 + 3r 4 g 2 + 5r 4 g| - 5r 8 ) 
Xl2 ~ 3f 8 - f 4 q\ - f 4 q\ - q\q\ ' 1 j 

X22 can be obtained by interchanging the indices of the expression of xn- Thermodynamic 
instabilities appear for all the black holes unless q\ — q 2 - The critical surfaces are listed 
below. 

3f 8 _ f 4 g2 _ ^4-2 _ ^2-2 _ q £ or q^-£ q 2 -£ 

3f 4 - q\ = 0, for q 1 ^0,q 2 = 0. (2.40) 
Critical exponents turn out to be same as (2.19). 
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3 Mean Field Analysis 



In the previous section, we studied the thermodynamic properties of R-charged black 
holes in five, four and seven dimensions. We have seen that, for generic values of charges, 
these black holes undergo continuous phase transition characterized by certain values of 
critical exponents. In this section, our aim is to provide a mean field description of this 
black hole phase transition. More specifically, within the grand canonical ensemble, by 
identifying the horizon radius as an order parameter, we provide Bragg- Williams like 
effective potentials 1 that are expected to decribe instabilities associated with these black 
holes 2 . On shell, this potential reduces to the grand canonical potential. Furthermore, we 
will see, near the critical surfaces, the effective potentials leads to the right set of critical 
exponents. 

This section has three subsections. In the first subsection, we analyze black holes 
in five dimensions. In particular, we provide the effective potentials for black holes with 
one or two charges. For the case of three non-zero charges, owing to the complexity of 
the thermodynamic variables, we have not been able to construct mean field potential 
in a simple closed form. In the next subsection, we carry out a similar exercise for 
black holes in four and seven dimensions. Finally, in the last part of this section, by 
trading horizon radius with charge density of the boundary gauge theory, we propose an 
effective potential which is expected to decribe the gauge theory in the presence of non- 
zero chemical potential. We illustrate this explicitly for simplest R-charged black hole, 
namely, for a singly charged black hole in five dimensions. 

3.1 Black holes in five dimensions 

We start with the simplest of R-charged black holes - one with a single non-zero charge in 
five dimensions. Thermodynamic quantities can be found out simply setting q~2 = qz = 
in (2.8, 2.9, 2.10). We suggest that it is described, by an effective potential 3 



Note that, as in grand canonical ensemble, T, [i\ are to be treated as external parameters. 
T also depends upon f, the dimensionless horizon radius, f will be treated here as the 

1 Similar constructions were also found useful in oreder to study various black holes and gauge theories 
around the Hawking-Page points. See for example [16-18]. 

2 Though in this paper, we will interchangeably use the words Bragg- Williams potential and mean-field 
potential, it should be recognized that BW potential and Landau mean-field potential are not the same. 
Close to T c , the later appears, in some cases, as an expansion of the former in terms of order paramer. 
We will comment on such a phenomenon in section 3.3. 

3 Constructional procedure is similar to the one described in [10]. 
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order parameter. As we will discuss below, our choice leads to the correct set of critical 
exponents near the critical line. We also point out that in JF, f can take any value; the 
relation between f and the mass and the charge of the black hole appears as a consequence 
of on-shell condition. 

The equilibrium condition of the system is represented by 

f = a (3,) 

This leads to 

f = (Af 2 - Jj 2 ) = q l + 2f 2 
2y/2iry/(2r 2 -p,l) 2n^/f 2 TYi 
This is the temperature of a singly charged black hole as can be seen from (2.8) after 
substitutions q\ = qs = 0. Now using (3.3) in (3.1), we get 

f_ f6 _ r 2 (g!+f 2 ) _ - 

87rGu 2 -2f 2 ) 16tt 1 ' 

So on-shell, T reduces to the negative of the pressure P which is same as the grand 
canonical potential. The extrema of (3.1) appear at 



= 0, 



\\JlA + 2vr 2 T 2 - 2^7T 4 f 4 - 7r 2 /x 2 f 2 , 

]- Jfil + 2n 2 T 2 + 2^/tt 4 T 4 - ir^T 2 . (3.5) 



2 

We note that within the range of our interest, namely for q x < 2f 2 (which leads to 
T > T c ) , fourth root of (3.5) corresponds to the only minimum. This equilibrium value 
of f increases with the temperature. 

Behaviour of T as a function of the order parameter for fixed Jjl\ and for different 
temperatures is shown in the figure (1). Clearly, the order parametr f changes continu- 
ously around the critical temperature. This is typical of a second order phase transition. 



The mean field potential T reporduces correct set of critical exponents discussed in 
sec. 2.1. Approaching the critical line along constant Ji\ = /2 lc , we see from the last root 
of (3.5) 

f - r c ~ (f - T c )5, with f c = -y^ic- (3.6) 
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-4 - 



Figure 1: T as a function of order parameter f for fixed jSj. The solid line is for the 
T = T c . The other two lines are for T > T c . As T increases, minimum of JF appears for 
larger values of f. From T = T c , the order parameter changes continuously with T. 



This leads to the critical exponent (3 = |. Similarly, the susceptibility, that follows from 
(3.5), behaves as 

x = ^L| f ~(r-r c )-i, (3.7) 

near the critical temperature. Therefore 7=3- The specific heat follows from (3.1). 

cj h = -r^| ffl ~(T-r c )-i, (3.8) 

which gives a = |. Finally, approaching the critical line with T = T c and using (2.10), 
we get from the last root of (3.5) 

f - f c ~ - (/2ic — Mi)^, with f c = ^-^^c- (3.9) 

where f c is the critical value of f at /2i = //i c . Hence, we have 5 = 2. We therefore 
conclude that the effective potential of (3.1) reproduces all the critical exponents given in 
(2.19). 

Next, we consider black holes with two non-zero charges. Thermodynamic variables 
can be found from (2.8,2.9,2.10) by substituting ^3 = 0. The mean field potential can be 
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expressed in the following manner: 

f: ( r , inlnl - 4r 4 ) r (12f 4 l/l?^ 2 - W 2 + /'T/'5) /Ttt 



Svr L 8(/i 2 -2f 2 )(/i 2 -2f 2 ) v /(/i 2 -2f2)(^_2f2 / 

(3.10) 

The potential is a function of the order parameter f and also depends on the external 
parameters T, Jx\ and // 2 . Note that in the limit Jl-i = 0, T reduces to that of (3.1). 

The equilibrium configuration is found by minimizing the potential with respect to 
f. This leads to 

f _ r(4f 2 -p.l-p%) r qi + q 2 + 2f 2 

27rV(^-2f 2 )(/i 2 -2f 2 ) 27r v /( f 2 + g 1 )(r 2 + g 2 )' 

This is the equilibrium temperature of the doubly charged black hole (2.8). Substituting 
(3.11) in (3.10), we get 

f = - 64,(^)5-^) = -isf<* + * )(f2 + « = (3 ' 12 » 

The values of f at which has extrema can be evaluated analytically. However, 
expressions are very large and unilluminating. We do not display them here. It turns out 
that in the region of our interest namely for 2f 2 — q\— g 2 > 0, there is only a single minimum 
of T . Typical behaviour of JF as a function of f for fixed pi and /2 2 is shown in the figure 
(2). We see that the order parameter f, at which the minimum appears, continuously 
increases as we increase the temperature. This shows that the phase transition is of 
second order. 

Various critical exponents can be evaluated as in the single charged case. Approach- 
ing the critical surface along constant fii and /i 2 , and expanding the minimum value of f 
close to the critical temperature, we get 

f-f c ~ (f -f c )i (3.13) 
Similarly, on constant T, p, 2 and T, Jii surfaces, we get 

^=wr {T - %rK (3 ' 14) 

respectively. The specific heat follows from taking second derivative of T with respect to 
T for fixed fii and Jjl 2 - This leads to 

C, ufi2 = -f~~(f-f c )-l (3.15) 

Finally, on constant T, // 2 and T, Jl\ surfaces, expanding around the minimum of T for fii 
and p, 2 close to their critical values, we get 

f-f c (fil-p>lc)*, f-f c ~ - (p 2 -p2c)K (3-16) 
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Figure 2: T as a function of order parameter f for fixed fti and /i2- The solid line is for 
the T = T c . The other two lines are for T > T c . As T increases, minimum of T appears 
for larger values of f. From T = T c , the order parameter changes continuously with T 

respectively. We now note that these lead to the same critical exponents that we got from 
our analysis in the previous section, see (2.19). 

3.2 Black holes in four and seven dimensions 

Here we reconsider the black holes in four and seven dimensions. As in previous subsection, 
our aim will be to understand their behaviours near the critical surfaces in terms of 
effective potentials. Though for all non-zero charges, it is in principle possible to construct 
the function T, it becomes computationally difficult to express T in terms of temperature 
and all non-zero chemical potentials in a closed form. We will therefore consider, in this 
subsection, black holes with a single charge. 

Thermodynamic variables for four dimensional holes with single charge can be found 
out by substituing g 2 = g 3 = g 4 = in (2.25). The effective potential can be written down 

as 




(3.17) 
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The saddle point of T gives the on-shell temperature of the black hole 

T= ^-ti = (3r- + 2g?)Vf 
4y/2iry/2r 2 - j% Au^f + ql 

This is same as what we would have gotten from (2.25) after setting three charges to zero. 

The critical line (2.31), while expressed in terms of T and pi reads as T c = y / 3/ii/27r. 
We will be interested in the rigion T > T c or equvalently in the region qf/r < 3/2. There 
is only one minimum of JF in this region and is given by 



f + ^ + ^4^-3^2. (3.19) 

Studying T around this point, as in the last subsection, we recover the critical exponents 
(2.19). 

For seven dimensional black holes, one can similarly construct the effective potential 
T . Following is the expression of T for single charge case. 



^6 



5r z -2n\ f- - I 1 



2v / 27rTi 



2f 2 -n\ v 2f2 - \A 



(3.20) 



As before at the saddle point this function reduces to grand canonical potential with 
temperature given by the one in (2.37) with q 2 = 0. The critical surface appears at 

T = \J\^-- For a given pi, this equation defines f c . In the range of our interests, namely 
for T > T c , there is only one minimum of T . This is given by 



f = V f + ^f- + \\j^f± - Wfif*. (3.21) 

Studying the behaviour of T around this point, we recover all the exponents that are 
given in (2.19). 



3.3 Proposal for gauge theory effective potentials 

We now wish to use the results of the previous subsections to propose an effective potential 
describing the gauge theory. While it might be possible to construct such a potential in 
complete generality, in this subsection we study the simplest case. We consider the four 
dimensioal finite temperature gauge theory where only one chemical potential is turned 
on. 

The right order parameter to work with, at the boundary, is the R-charge density 
p\. The conjugate chemical potential is jii. We can now use the first equation of (2.10) 
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to express (3.1) in terms of p\. This leads to the following equation for T: 




8tt 



s]a^pI-N^\ Pi . (3.22) 



In writing this down, we have reinstated G the five dimensioal gravitational constant, 
various factors of I and also used the relation G = 7r/ 3 /2iV 2 . Here N c is the number of 
colours. 

One can now check that, for T > T c = 7r//i lc , the minimum of T occurs at 



= (NW\ {/if + 2( 7 r 2 r 2 + y 7 r^-^ 7 r^)} 2 
Pl V 32tt 2 J -tf + 2(vr 2 T 2 + v^T 4 - p 2 7r 2 T 2 ) ' 1 ' ' 

Now expanding (3.22) near T = T c , for fixed we get 

Pl - Plc ~(T-T c )5, (3.24) 

where the critical charge density is given by p lc = 9pl c N^ /(327r 2 ). The above equation 
leads to P — 1/2. Similarly, one can find the other critical exponents given in (2.19). 

To this end, we would like to make a couple of comments. Firstly, one may wish to 
expand T near T = T c in powers of (p — p c )/p c t° get a Ginzburg-Landau potential as 
an expansion in terms of order parameter. However, in doing so one does not reproduce 
the right exponents. Higher and higher powers seem to conspire in order to recover 
the results of (2.19). Secondly, since the gauge theory is a strongly coupled, directly 
computing the potential (3.22) seems a difficult task. However, in [2], an attempt was 
made to construct Ginzburg-Landau potential for the system via a regulated field theory 
model. Unfortunately, all the critical exponents predicted by supergravity configuration 
could not be recovered. 

We end this section with a note that, at least for single charge holes, similar gauge 
potentials can be easily constructed in three and six dimensions simply using (3.17) and 
(3.20) along with the thermodynamic relations listed in the previous section. 



4 R-charge Conductivity 

So far we have discussed equilibrium thermodynamics of the black holes and their dual 
gauge theories. In this section we will turn towards non-equilibrium thermodynamics of 
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this system; in particular we will compute the R-charge diffusion constant using gauge- 
gravity duality. There are two ways to compute R-charge diffusion constant. One can 
consider retarded correlation function of R-current and look for a pole in the g 2 -plane 
where q represents the momentum. Or one can make use of Green- Kubo formula to 
calculate the R-charge conductivity A and find R-charge diffusion coefficient D from the 
relation A = D . x, where x is the susceptibility matrix. If the susceptibility matrix is 
already known, the advantage with the the second method is one need to compute the 
correlator only at the zero momentum limit and upto first order in frequency. Since, in 
the present text we have already evaluated the susceptibility matrix for various cases we 
will consider the latter method only. 

The procedure of computing R-charge diffusion constant using Kubo formula is 
explained in detail in [19,20]. For convenience of the reader we recount the essential steps 
involved in the following. For this method we need to find out the retarded correlator of 
the R-current: 

oo 

&£>(u,$) = -ij dte^9(t) J dxe-** ([J x (t,x),J x (0,0)]) • (4-1) 

— oo 

at the zero momentum limit. The R-charge conductivity A is given in terms of retarded 
correlator by Green-Kubo formula: 

A= lim ^W = o)] 

oo 

= lim^- J dte tuJt J dx ([J x (t,x) : ,7.(0,0)]) . 

— oo 

The retarded correlator can be computed as follows: We find perturbations in various 
modes around the black hole solution that satisfy linearized equation following from equa- 
tions of motion and the 'incoming wave boundary condition' at the horizon. Then we 
calculate the boundary action which can be obtained by plugging in the solutions in the 
expression of total action evaluated at boundary. The total action consists of three pieces: 
the bulk action, the Gibbons-Hawking term and the counterterm required to cancel the 
divergences. The retarded correlator can be obtained by taking derivative of the boundary 
action with respect to boundary value (of the perturbation in the gauge field) twice. 

Let us demonstarte it in the case of five dimensional black hole with three charges. 
For the present purpose it is sufficient to consider perturbations in the tensor (metric) 
and the vector (gauge fields) modes around the black hole solution and keep the scalars 
unperturbed. So perturbations are of the form: 

<fo* = gjg + V » 4 = Aj°>+4, < = 1,2 . (4.2) 
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In order to determine R-charge conductivity it is enough to consider perturbations in (tx) 
and (xx) component of the metric tensor and x component of the gauge fields. Moreover 
one can choose the perturbations to depend on radial coordinate r, time t and one of the 
spatial worldvolume coordinate z. 

A convenient ansatz with the above restrictions in mind is 4 : 
h tx = g 0xx T(u) e~^\ h zx = g 0xx Z(u)e-^ t+ ^, A\ = | A\u) e~^ z . (4.3) 

Here u and q represent the freequency and momentum in z direction respectively and we 
set perturbations in the other components to be equal to zero. 

Next step is to find linearized equations following from the equations of motion at 
zero momentum limit that the perturbations has to satisfy. It turns out that in linearized 
equation at zero momentum limit metric perturbation Z(u) decouple from the rest. One 
can further eliminate T(u) reducing it to equation for perturbations in gauge field only. 
For five dimensional black hole with all the three charges the equation becomes in this 
notation 



//' H[ H' 2 H' 3 \ ., u 2 U Al u(l + h)v /n , W1 i \ a 



+ fc 2 (l + h)(l + k 3 )A 2 + fc 3 (l + + k 2 )A 3 



(4.4) 



0, 



while equation for the other two components can be obtained by cyclically permuting 
(123). 

We require solutions to these equations with "incoming wave boundary condition" 
at the horizon at u — 1. Since f(u) vanishes at u — 1 one can solve the indicial equation 
near u — 1 and obtain two solutions which have asymptotic behaviour near u — 1 as 
(u— \) v with v = ±itu(T /2T). The incoming wave function corresponds to the negative 
value and so incorporating the boundary condition corresponds to the following ansatz: 

f-iui(T /2T) 

A i = J —— ai(u), i = l,2,3. (4.5) 

Since in the sequel we need solutions upto first order of u only we can write it as 

di(u) = Pi(u) + iusi(u) + 0(co 2 ). (4.6) 



4 Here is a comment about notation in this section. To facillitate comparison with expressions in 
existing literature we have made a few changes in our notation. In what follows, instead of r we use 
u = (r + /r) 2 so the boundary is at u — and horizon is at u = 1. Similarly the charge parameters qt 
will be traded for hi = {qi/r\) and we will also introduce To = ('"+/ 71 ")- In this notation the solution for 

m 

black hole is given by Hi = (1 + kiu), H — TT Hi where m is number of charges. For five dimensional 

i=i 

m m 

black hole / = H — T7 (1 + ki)u 2 and for four or sweven dimensional black hole / = H — pj(l + ^i) u3 - 

i=l i=l 

We will also use Co = oj/(2ttT ) 
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In terms of a^'s the boundary action appears to be 

N 2 T 2 f f 

Sboundary = lim c / dtdf — [H 2 a[ai + H 2 a' 2 a 2 + H 2 a 3 a 3 ], (4.7) 
m^o lb y /t 

where we have kept terms quadratic in perturbation. The retarded correlator can be 
obtained by taking derivative of the total action with respect to boundary value twice. 

Proceeding further with this generic charge leads to expressions, which are techni- 
cally complicated. Moreover, one unusual problem here is that we cannot obtain special 
cases, such as two charge solution, starting from generic cases (three charge) by setting 
one of the charges to zero. This happens because obtaining these equations involves rescal- 
ing of fields A l (u) with respective charges and forces us to do the analysis in each case 
separately. Therefore, we will restrict ourselves to the R-charge conductivity (for black 
holes in five, four and seven dimensions) in the simplest cases only, deferring analysis of 
more complicated cases to the appendix. We will also adjoin a brief outline of the general 
method that has been used to solve all these cases. 



4.1 Five dimensional black hole with two charges 

In this case we have two U{\) gauge fields. The solution is 

#i = (l + M, H 2 = (l + k 2 u), H = H X H 2 f(u) =H-(l + k 1 )(l + k 2 )u 2 . (4.8) 
The perturbed equations are given by, 

A „ ( f W H[\ Al Co 2 H A 



(4.9) 



fH'i 



2 



-(fci(l + fc 2 )Ai + fc 2 (l + h)A 2 ) = 0, 



'f_H' , n H^\ A/ , Cj 2 H 
u(l + k 2 



A 2 + \ — + 2— ) AL + — f A 2 

- \f n h 2 J 2 up 



(4.10) 



r (k 1 (l + k 2 )A 1 + k 2 (l + k 1 )A 2 )=0. 
3 tii 



l 2 

. We observe that the equation for A 2 can be obtained from equation for A\ by inter- 
changing ki,Hi,Ai with k 2 ,H 2 ,A 2 respectively. From now on for convenience, we will 
write perturbed equation for only one of the gauge fields. The ansatz is given by 

f-iui(T /2T) 

A i = J —— ai (u), i = l,2,3. (4.11) 

We need to solve for a^u) in terms of the boundary value of the gauge fields. Since in 
the sequel we will need solutions upto first order of u> only we can write it as 

fli(w) = [Pi(u) + iwsi(u) + 0(Co 2 )] . (4.12) 
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where % takes value upto 2 for this particular case. It turns out from (4.9), that p\ and si 
satisfy the following equations respectively: 



KH + (y-^) P'i( u ) 

h (f H'\ . . 1 + fci 



h{l + k 2 ) k 2 (l + h) 

-Pi H 77 V2 



Hi H 2 



o, 



fcl /7' W'\ , v 1 + fcl 



fci(l + A;2) , fc 2 (l + ^i) 

■Sl H 77 s 2 



0. 



Similarly by interchanging k±, pi, Si with k 2 , p 2 , s 2 we obtain equation for p 2 and s 2 . 
According to our boundary condition all these solutions need to be regular at u — 1. 
Without any loss of generality we can choose = . 

The solution for Pi(tt) are given below: 



/ \ , . faUv , hl + k 2 

P2(u)=b W (l + —)-b 2Q -—U, 



(4.13) 



where foio and 620 are boundary value of A\ and A 2 respectively. The solutions of Si(u) 
are more involved and given by 



si(u) = hiu + (h 2 + h 3 u) log [1 + (1 + hi + k 2 )u], 
s 2 (u) = hu + (l 2 + l 3 u)log [1 + (1 + h + k 2 )u\. 



(4.14) 



Note that at u = 0, Si(u) = as mentioned before. The constants are given by 

hi = b 10 Ji + 620^20, h 2 = b 10 L w + b 20 L 20 , ( 4 -15) 



where 



J\o — ki 
J20 — ~k 2 



{k\ - k\{k 2 - 3) + k 2 {k\ - k 2 - 6) + k 2 {k\ + 3A; 2 + 2) 

4(1 + h)(i + k 2 ) (i + h + k 2 ) 

k\ - k\{k 2 - 3) - h(k 2 2 + 2k 2 - 2) + k 2 {k\ + 3A; 2 + 2) 



The values of L w and L 20 are 

-£ao - 
L 2 o = 



4(l + A; 2 ) 2 (l + A; 1 + fc 2 ) 



2 + 2kj + 2k 2 + k\ + fci(4 + 3A; 2 ) 
(1 + fci + A; 2 ) 2 
(1 + fci)fc 2 (2 + fci + fc 2 ) 
(l + fc 2 )(l + A;i + A; 2 ) 2 



(4.16) 
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h 3 is given by 



fc 2 (l + fci) ^ 

1 + fc 2 



/i, / 2 and / 3 are given by 

^1 = ^10^10 + ^20^20, ^2 = ^10-^10 + ^20-^20, 



(4.17) 
(4.18) 



where Jio(^2o) can be obtained by interchanging fci with fc 2 in J 2 o(</io) and Lio(£2o) can 
be obtained by interchanging k\ with fc 2 in L 20 (-fcio) 

The second step consists of evaluation of the boundary action quadratic in the 
boundary values. The boundary action turns out to be of the following form. Since we 
need only imaginary contribution of 0(a)) we keep the other terms implicit. 



'boundary 



lim 



u^o 16 
=0(u> ) + 



dtdx — [H 2 a[ai + -f/f a 2 a 2 ] 



N 2 C T \ 
16 



-iuo 



'10 



Jio + Jio(l + h + fc 2 ) - (T /2T)(fci + fc 2 ) 



+b w h 



10^20 



^20^20 + (1 + h + fc 2 )(J 20 + L 20 ) 



+b 



20 



L w + L 10 (l + fci + k 2 ) - (T /2T)(fc! + fc 2 ) 



+ 0(c^) 



From the above boundary action we obtain the following components of the retarded 

J L0 + Jio(l + h + k 2 ) - (T /2T)(h + fc 2 )] + G(cu 2 ), 

^20^20 + (1 + fcl + k 2 )(J 20 + ho) + 0(U 2 ), 

Lio + Li (l + fci + fc 2 ) - (T /2T)(fc! + fc 2 ) 



Green's function 






GSf ) (cu,g = 0) = 


-0(u ) 


■ ~N 2 T 2 \ 
lU 8 [ 


GSf } (^g = 0) = 


-0(Q°) 


■ ~N 2 C T 2 r 
- W 16 I 


G£ R) (u,q = 0) = 


=0(a>°) 


.~^t 2 r 
w 8 [ 



o(^ 2 ). 

(4.19) 



The components of the R-charge conductivity then obtained as 



An =2. 



N'T 2 



16 

Al2 -^6~ 
N 2 T 2 



A 22 —2. 
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Jio + Jio(l + fci + fc 2 ) - (To/2T)(fcx + fc 2 ) 
^20-^20 + (1 + fci + k 2 )(J 20 + L 20 ) 
Lio + L W {1 + fci + fc 2 ) - (T /2T)(h + fc 2 ) I . 



(4.20) 



The susceptibility matrix Xi? can be written in this notation as 



Xij 
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k\ + fcifc 2 + fc 2 - 5fci - 2 



v / fc7fc 2 "(fci + fc 2 ) 



Vhfa(h + fc 2 ) fcf + fc 2 fci + fci - 5fc 2 - 2 y fci + fc 2 - 2 



(4.21) 
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The components of the diffusion coefficient can be obtained by using the relation T> = 
X-X 1 - Since all the components of susceptibility diverges on the critical line k\ + k 2 = 2 
and conductivity remains finite it follows that all the components of the diffusion coeficient 
vanish. Vanishing of the diffusion coefficients is consistent with slowing down at critical 
point. 



4.2 Four dimensional black holes 

Next we consider the black holes in three dimensions with two charges (single charged 
case has already been discussed in [7]). There can be maximum four charges as there are 
four commuting f/(l)s. Leaving the more general cases for the appendix, here we consider 
black holes with two £/(l)s. The equation for gauge field is given by 

f W H[\ Al Co 2 H A 

T - W7 + V M + —el A * 
f n Hj uP 

v 2 C\+kA 1 ' 

~ f rr2 W + k2)A, + fc 2 (l + h)A 2 ) = 0. 

Equation for A 2 can be obtained by inchanging (Ai, k\,H\) with (A 2 , k 2 , H2) respectively. 
The ansatz consistent with the boundary condition is 

f-iu(T /2T) 

= i a i( u )i a i( u ) = Pi( u ) + iusi(u). (4.23) 

1 "T" K%Uj 

Using above ansatz we obtain equations for pi,p 2 , si, s 2 . The solution for pi is given by, 

/ % . 2fci k 2 (l + h) 

pi(«) = 6io(l + —u) + b 20 3( \ + fc2) V (4.24) 

while p 2 can be obtained from above expression by permuting k\ and k 2 . The equations 
for Si are more involved and can be written as 



#1 



2 r 



h(l + k 2 ) k 2 (l + h) ' 



(4.25) 



- 2vp 1 (u)'j + vpxiu) (Jj^ ~ y ) = 0, 



where we used v = T /(2T). Equation for s 2 can be obtained by permutation, as usual. 
Solution is given by, 



si(u) =h\u + (h 2 + h^u) arctan 

/ 

+(h 4 + h 5 u) log(- ), 

1 — u 



^/3-k 2 + 2k 2 -k 2 + 2fci(l + k 2 )u 
2 + (1 + £4 + k 2 )u 



(4.26) 
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and similarly for s 2 (u). The parameters are as follows: 

h = ( (3 + 2*4 + 2k 2 + 2k x k 2 ) \ - 2(1 + A^A^ + A; 2 (l - k 2 ) - hk 2 (l + k 2 ) + k 3 2 (l + k 2 )) 



>2Q 



+ fcx(l + fc 2 ) Ak\ - kj(7k 2 - 4) + (fc 2 + l)k 2 2 + hk 2 (2k 2 - 7) 

v 



no 



x 



&2 = — 



9(1 + fc!)(l + k 2 ) 2 [-3 + k 2 - 2k 2 + k 2 2 - 2k x - 2k x k 2 \ ' 
2h 3 , (1 + A;i) 



fci fci(l + fc 2 ) 

/i 3 = (2A;i(l + A;2)[9 + ^ + +3A;2 + ^ + ^ + ^(5A;2 + 13) + A;i(21 + 8A;2-A; : 



(4.27) 
(4.28) 



fc 2 (l + fci) 9 + fc? + 21fc 2 + 13A;| + k\ + A; 2 (13 + 5A; 2 ) + fci(21 + 20A; 2 + 5A: : 



>10 



'20 



3(1 + fc 2 )(3 - A; 2 + 2A; 2 -k% + 2k x + 2k l k 2 f/ 2 ' 
/i 4 = (3/2)u6io, 

/i 5 = kivb w - A; 2 + , X > vb 1Q , 
2(1 + fc 2 ) 

where we have used v = (T /2T) and the parameter / 3 : 



(4.29) 
(4.30) 
(4.31) 



h = (2(1 + h)k 2 [9 + k\ - k\{k 2 - 1) + 21A; 2 + 13k 2 2 + k\ + fci (3 + 8A; 2 + 5A; 2 )] b 20 



fci(l + k 2 ) [9 + kf + 21k 2 + 13A; 2 + k% + kf(13 + 5k 2 ) + 21 A* + 20£; 2 + 5fc|] 6 
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3(1 + fc0(3 - A; 2 + 2A; 2 - fc| + 2^(1 + A; 2 ) 3 / 2 
The boundary action is given by 



(4.32) 



= lim 



N~ 2 T 



J dtdx [H 2 a[ai + i/f a 2 a 2 ] . (4.33) 



Oboundary ^ ^ 

Now we can follow the same procedure as in the d — 5 case and obtain the R-charge 

conductivity by differentiating twice with respect to 610 and 6 2 o- The (11) and (12) 
components of R-charge conductivity turn out to be 



A 



11 — 



N£ T Q v 2 [4kf{2 + k 2 ) + 9(3 + 3A; 2 + k 2 ) + fc?(45 + 3lk 2 + k 2 ) + h(63 + 54k 2 + llfcf)] 



A 



3QV2 
n! T v 



9(l + fci)(l + A; 2 ) 
k%(2 + k 2 ) + A; 2 (9 + 9A; 2 + 2A; 3 ) + kf(2 + 18k 2 + \3k\ + 2k\) + 



36v/2 9(1 + fc!) 2 (l + fc 2 ) 2 
k 2 (9 + 54k 2 + 48A; 2 + 13fc| + k 2 ) + h(9 + 54k 2 + 54A; 2 + 18A;| + 2fc|)] . 



(4.34) 
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The (22) component can be ontained by interchanging k±, b w with fc 2 , 620 respectively. 
Once again since all the components of susceptibility diverges near critical point, while 
conductivity remains finite, as evident from the above equations, the diffusion coefficients 
will vanish. 



4.3 Seven dimensional black hole with two charges 

The seven dimensional case corresponds to the M5 brane. Since the transverse space is 
five dimensional only two U(l) charges are possible. The equations for the perturbation 
in gauge field is given by 

A " +(t-^ + 2 ^--) A> i- + <u2 U * [Ml + fcMi + Ml + k i)M = 0. (4.35) 

\J ft tlx UJ J tlx 1 J 

The other equation can be obtained by interchanging Ax, kx with A 2 and hi respectively. 
The ansatz appropriate to the boundary condition is 

t — iu>v 

Ax(u) = J ——(px(u)+iusx(u)). (4.36) 
tlx 

The explicit form of the solution to Sj(w)'s are obtained but they are not written here 
as they are quite complicated. But the conclusion remains the same i.e. as we approach 
critical line diffusion constant goes to zero. We explain in the appendix (see last section) 
why two-charge case for M5 brane is complicated than the same for D3 and M2 brane. 



5 Conclusion 



In this paper we have investigated some of the features of equilibrium and non-equilibrium 
thermodynamic properties of R-charged AdS black hole in its full generality near the crit- 
ical point. Furthermore, by making use of the gauge-gravity correspondence, we tried to 
extract some properties of strongly coupled boundary gauge theories at finite temperature. 
More specifically, we analyzed five, four and seven dimensional black holes with multiple 
charges. This led to the construction of Bragg- Williams potentials associated with these 
black holes. We discussed how Bragg- Williams potential allow us to make a proposal for 
the gauge theory effective potential (above the critical temperature). This effective po- 
tential is usefull in determining the equilibrium properties of the strongly coupled gauge 
theories in the presence of non-zero chemical potential. Subsequently, we analysed some 
of the non-equilibrium properties of the black holes and its duals. Though, for single 
R-charged black holes, thermal properties were discussed in [1] and [7], we extended these 
results in several ways. Besides providing a general framework to solve the perturbed 
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gravity equations, we determined the conductivities and diffusion coefficients for multiply 
charged black holes. 

Our work may be extended futher by considereing rotating R-charged black holes. 
They were consuructed, for example, in [21]. Another question that is of crucial impor- 
tance is to understand the fate of these R-charged black holes below the critical tempera- 
ture. As we saw, these holes become unstable. However, the new phase to which it crosses 
over is, to our knowledge, not yet clearly understood. 



Acknowledgements 

We have benefitted from useful conversations and communications with Somen Bhat- 
tacharjee, Sankhadeep Chakrabarty, Sumit Das, Alok Kumar, Kengo Maeda, Makoto 
Natsuume, Takashi Okamura, Binata Panda, Kalyana Rama, Bala Sathiapalan, Sreeku- 
mar Sengupta, and Goutam Tripathy. We thank organizers and participants of ISM 08 
at Pondicherry where part of this work was presented. 



A General solution to linearized equation for pertur- 
bation 

In this appendix we briefly sketch the structure of the solutions and the general method 
that have been used to solve linearized equations for the perturbations in the gauge field. 
The equations satisfied by Pi and Sj in each of the cases, i.e. for five four and seven 
dimensional black hole can be written in the following general form: 



(A.l) 



(A.2) 

where m is the number of R-charges. The parameters c, n, b and m have following values 
for the different cases: 
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• Five dimensional black hole: c = 0, n — 1, b — 1, m — 1, 2, 3. 

• Four dimensional black hole: c = 0, n = 2, b = 1, m = 1, 2, 3, 4. 

• Seven dimensional black hole: c = —1, n — 3, b — 4, m = 1, 2. 

The solutions should be regular at horizon ie. at -u = 1. As already mentioned in the 
main text, without any loss of generality we can also impose the condition that at u — 
the term linear in frequency should vanish. 

Solution for p^s can be written generically as 

Pi(u) = b i0 + do u l , (A. 3) 

where Qo depends on all the k^s and b^s. Clearly b^ is the constant representing bound- 
ary value of Ai(u). For five and four dimensional black hole I = 1 while for seven dimen- 
sional black hole 1 = 2. 

Now we turn to solution for Sj's. The function f(u) in the metric generically can be 
written as, 

f = (l-u)F(u), (A.4) 

where F(u) is a polynomial in u of order v, and so has v number of zeroes given by 
u = Uj,j = 1, 2, .., v. For four and five dimensional black hole v — m — 1 while for seven 
dimensional blackhole v = m + 2 

Then most general solution for Sj can be written as 

V 

Si = Pi{u) + Qij(u) log(l - «/uj), ( A - 5 ) 
j'=i 

where Pi(u) and Qij(u) are two sets of functions. In order to ensure Sj(0) = 0, we have 
-Pj(O) = 0. We require our solution to be regular at the horizon so log(l — u) term is 
absent. This is the genral form representing the solution for various cases as discussed 
in main text as well as in appendix. For each particular case we need to determine the 
functions Pi(u) and Qij{u). 



B Five dimensional black hole with three charges 

In the case of five dimensional black hole with three charges the values of the parameters 
are c = 0, n — 1, 6 = 1, m = 3 and I = 1. Since / = (1— u)(l + (l+ki+k2+k 3 )u— k^k^u 2 ) 
we get 

F(u) = (1 + (1 + h + k 2 + k 3 )u - k^hu 2 ), (B.l) 
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which is a quadratic polynomial implying v — 2. It is helpful for calculation if we write 
the second term in (A. 5) as 

Qi2 



Qalog(l-u/ Ul )+Q t2 \og(l-u/u 2 ) = Q * + Qa i og (F( M )) + ^ 



Using I — 1 we write down from (A. 3) as 

= b i0 + c i0 u. 

Cio turns out to be 



log 



/ M 2 («l ~ It) 

\Ui(u 2 - u) 

(B.2) 



(B.3) 



ClO 



fci&io - (1 + h 



ho + 



'30 



1 + A; 2 1 + A; 3 

Other components of Qo can be obtained by cyclic permutation of (123). 
Expressions of Sj(-u) are as follows: 

f \ u \ i u \ u w f^ + ^ + h + k 2 + k 3 + ^A)u 

si(u) = h x u + (/la + h 3 u) log — — — - ? =— 

\2 + (1 + fci + fc 2 + % - V A)u 

+ (/i4 + hu) log[l + (1 + fci + £; 2 + k 3 )u - kik 2 k 3 u 2 ], 



(B.4) 



while 



A 



(1 + h + k 2 + k 3 f + Ak x k 2 k 3 . 



(B.5) 



(B.6) 



hi for i = 1,2,3,4,5 are constants, which depend on k±, k 2 and k 3 and boundary values 
of the gauge fields, namely, b w , b 20 , b 30 . Introducing some additional parameters, hi can 
be written succintly as 



hi 



(2 + ki + k 2 + k 3 - hk 2 k 3 )v 
4(l + fc 2 )(l + A; 3 )A 



hJi 



ho - 



k 2 J 2 



ho ~ 



k 3 J 3 



^30 



(B.7) 



1 + h ±u 1 + A; 2 l + k 3 
These new parameters Ji, J 2 and J 3 are fourth order polynomials in fcj's and are given as 
Ji - fc?(-2 + A; 2 + fc 3 ) - fc*(A; 2 - k 3 f + 4(A; 2 + fc 3 )(^ + k 2 k 3 ) +k\- 4k 2 k 3 

J 2 =[fc3 + ^ _ ^,3 _ ^,2^ _ ^2 + 3fc 2^ 3 + %k 2 kl + 10£;i£; 2 £; 3 ] 
+ [2(-fc? + k\ + 6A:|) + 4(fci + fc 2 )/c3] + [fci + k 2 + 7fc 3 ], 

J 3 =[fc3 _ ^3 + fc 3 _ 3£2£ 2 + g^2 + _ 3^2 _ fc 2^ _ ^2 + lOA^Afefcg] 

+ [2k\ + 6^ + 2A| + 4^3 + ±k x k 2 \ + [k 3 + 7k 2 + h}. 
Similarly h 2 is given by, 



(B.8) 



h 7 



1 + fci 



2v^ 

+ (l + A; 2 )(l + A; 3 )L 1 & 1 o 



2fci(l + fc 2 )(2 + fci + k 2 + fc 3 - hk 2 k 3 ) 



koL 



2-^2 



1 + k 2 



} 20 



k 3 L 3 
l + k 3 



b. 
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(B.9) 
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where Lj's are 

Li = -1 + k\ + fa + 3kj + kl + 3k\ + 3k^k 3 — kf + 2k\k 2 k 3 + k 2 + Qfafa + 3k 2 k\ - 3fa + fak\ 

+ V6fak 2 k 3 + <6k\k 2 k\ + fak\ + §fak\k 3 + k\[-3 - k 3 + 2k 2 2 k 3 - fa + 8k 2 k 3 + 2k 2 k% 
L 2 = (l + fa + k 2 -k 3 ), 
L 3 = (l + fa-fa + fa). 



h 3 is given by 



3 2 2 2(1 + A; 2 )(1 + A; 3 )A 3 / 2 



{k 2 M + k 3 N), 



(B.10) 
(B.H) 

(B.12) 



where v depends on fa's in the following manner 

V(l + fci)(l + fc 2 )(l + fc 3 ) 
2 + fci + k 2 + £; 3 - ^1^2^:3 

M is given as follows. 

M — — 2k 2 (l + k 2 )(2 + fci + A; 2 + fc 3 - fcifeA*) M 1 + fc 3)(l + fa + k 2 - A; 3 )6i 
+ fc 3 (l + fcx)(l - fa + k 2 + k 3 )b 30 + 

(1 + fcx)(l + fc 3 ) [ - 1 + k\ - k\ + k 3 + 3k\ + A;3 - 3^ - fc 3 A| 

+ k\(3 + k 2 + 3k 3 + 6k 2 k 3 + 2k\fa) + fc 2 (-3 + fc|)+ 

fci(l + 6A; 3 + 2k\k 3 + 3^ + 16A; 2 A; 3 + §k 2 k\ -k\ + Sfakj + 2^^)J b 20 . 

N can be obtained by interchanging k 2 and k 3 in the above expression. The rest of the 
parameters and h$ are 



(B.13) 



hi = yOio 

h = ^-(l + fa) 



fa 



1 + fa 



'10 



l + k 2 



J 20 



l + fa 



A30 



(1 + fa) 



(B.14) 



Evaluating the boundary action 

^boundary 



N^T 



36^2 



lim J dtdx ^ (Hfa^ai + H 2 a 2 a' 2 + H 3 a 3 a f 3 ) , 



bon,b 3 



(B.15) 
This 



one obtains the expression of the boundary action upto quadratic in 6 10 , u 2 o,u 3 q. 
expression is too long and so we are not writing it explicitly. The R-charge conductivity 
can be obtained by differentiating boundary action twice and are given by 



An = 



N^ /2) T 
36^ 



k((-l + k 2 k 3 ) + k\(fa 2 k - 3 - 2(5 + 3A; 3 ) + fc 2 (-6 + fc|)) 



- 4(2 + 2A; 3 + kl + k 2 2 (l + fa) + A; 2 (2 + 2A; 3 + k 2 3 ) - fa(16 + Ufa + 5k 2 3 + k 2 2 (5 + 4fa) + 



2k 2 {7 + 5k 3 + 2k 2 3 



2(1 + fa)(l + k 2 )(l + k 3 y 



(B.16) 
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A 



12 



7V C (3/2) T 1 
36V2 4 

f 

(l + fci) 2 + 



fci(4 + h + k 2 + fc 3 )(-2 -k 2 -k 3 + h(-l + k 2 k 3 )(l + k 3 )) 



(1 + fc 2 ) 



(B.17) 



Other components of the conductivity matrix can be obtained from above two com- 
ponents by cyclically permuting k±, k 2 and k 3 . The diffusion coefficients can be obtained 

3 

from the relation Dy = ^2 XikXkli where \ is the specific heat. Since all the components 

k=i 

of specific heat diverges the diffusion coefficients turn out to be zero, as expected. 



C Four dimensional black hole with three charges 



For four dimensional black hole with three charges we have c = 0,n = 2,b = 1, m = 3 and 
/ = 1. From f(u) = {l-u)[l + {l + k l + k 2 + k 3 )u + {l + k l + k 2 + k 3 + kik 2 + k 2 k 3 + kik 3 )u 2 ) 
we get 

F{u) = 1 + (1 + h + k 2 + k 3 )u + (1 + h + k 2 + k 3 + k x k 2 + k 2 k 3 + k x k 3 )u 2 ), (C.l) 
which is quadratic. The solutions are given by 

(l + fci) 



2k\ 

Mu) = b m (l + — > + 3 (1 + fc)(1 + iW 



k 2 (l + k 3 )b 20 + k 3 (l + k 2 )b 30 



(C.2) 



while other p^s can be obtained from above expression by permuting (123) cyclically. 

Using (A. 5) Sj(u)'s is given by, 

si(w) = h x u + (/i 2 + h 3 u)A(u) + (/i 4 + /i 5 u) logF(-u), 
s 2 (u) = hu + (l 2 + l 3 u)A(u) + (U + hu) logF(u), 
*3(w) = 777.1^ + (m 2 + m 3 -u)A(w) + (m 4 + m 5 w) logF(w), 

\/3 - k\ - k\ + 2£; 3 - k\ + 2fc 2 (l + fc 3 ) + 2^(1 + k 2 + fc 3 ) 



(C.3) 



A(u) = arctan 



2 + (1 + h + k 2 + k 3 )u 



-u 



Here once again we have combined two logarithms as in (B.2). The values of the param- 
eters are as follows: 



h 3 
h 5 



k 3 - 2fci(l + k 2 ){\ + k 3 ) \h 2 + (1 + ki)k 2 {\ + A; 3 )/ 2 + (1 + k x )k 2 k 3 m 2 



3(l + fc 2 )(l + fc 3 ) 

2fci(l + k 2 )(\ + fc 3 )5 10 + (1 + k x )\k 2 {\ + k 3 )b 30 + fc 3 (l + k 2 )b 3Q ] 

2(l + k 2 )(l + k 3 ) 



(C.4) 
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h 2 = ^4(1 + fci) 3 + 2*4 + 2£; 2 + 2k 3 + 2k x k 2 + 2k 2 k 3 + 2k 3 h 

fc 2 (l + fc 3 )(l + fci + fc 2 - fc 3 )&20 + M 1 + fc 2 )(l + h-k 2 + k 3 )b 30 

- (1 + fc 2 )(l + fc 3 )[4 + # + + ki + 4A; 3 + 6k 2 3 - 2kl(l + k 3 ) + 2k\{7 + 3A; 2 + 3A; 3 ) + 
2A; 2 (-3 - k 3 + fcf ) - 2A; 2 (-3 - 3A; 3 + kj + fcf) + 2^(17 + + llfc 3 + A$ + llfc 2 + 3A; 2 A; 3 ) 

- 2fci(-15 + £; 2 - llfc 3 + A$ + A$ + k\ + klk 3 - llk 2 - Ak 2 k 3 + k 2 k\) 



x (1 + k 2 ){\ + k 3 ) k\ + kl- kj(l + k 3 ) + (k 3 - 3)(1 + k 3 f - k\{\ + k 2 + k 3 ) 



- k 2 (5 + 6A; 3 + fcjj) - fci(5 + fc 2 + 6£; 3 + fc 2 + Qk 2 + 6fc 2 £; 3 

x ^(3 - k\- k\- k1 + 2ki + 2k 2 + 2k 3 + 2k x k 2 + 2A; 2 A; 3 + 2k 3 k x ) 

(C.5) 

where we have used v = T /(2T). The components of R-charge conductivity turn out to 
be 



An = 2 [Ak\ 



2 + k 2 + k, 



x 



+ k\ 53 + hk 2 2 + 43A; 3 + 5k 2 3 + k 2 (A3 + 23k 3 ) 

+ k\ 108 + k\ + 130A; 3 + 43£| + fcf + ^(43 + 26A; 3 ) + 2A; 2 (65 + 60A; 3 + 13fc 2 ) 
+ 9 3 + 6A; 3 + Ak\ + k\ + fcf (1 + k 3 ) + fcf (4 + 6A; 3 + 3kj) + k 2 {6 + 10k 3 + 6kj + k 3 3 ) 
fci [90 + 144A; 3 + 74A; 2 + llfcf + fc|(ll + 10fc 3 ) 

+ A; 2 (74 + 82A; 3 + 29A; 2 ) + A; 2 (144 + 187£; 3 + 82k 2 3 + 10fc|) 



A 



9(1 + fci)(l + fc 2 )(l + fc 3 )(l + k l + k 2 + k 3 )J' 
k 2 + /c 2 /c 2 + fci(l + 4A; 2 + k 2 2 ) 



12 



9(1 + fci) 2 (l + A; 2 ) 2 (l + fc 3 )(l + h + k 2 + k 3 



2/co 



2 + h 



+ 2k\ 



(C.6) 



2 + k 2 + k 3 \ + 



k\ 



k\ 



22 + 14A; 3 + k 



- 2 



-9-9k 3 + ki + k 



- 36 - 26k 3 - Ak( + k 



+ 



22 + 4A; 2 + 14A; 3 + k 2 3 + k 2 (20 + 7A; 3 ) 



9-9k 3 + k 2 + fcf 



- 36 - 26fc 3 + 4fc 2 + ^ 



+ k\ 22 + Ak\ + 14fc 3 + k 2 3 + fc 2 (20 + 7fc 3 



+ fci 



36 + 2kl + 26fc 3 - 4fc 2 - k\ + fc 2 (20 + 7fc 3 ) + fc 2 (50 + 24A; 3 - 8kj 



(C.7) 



while other components can be obtained by cyclically permuting ki, k 2 and k 3 
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D Four dimensional black hole with four charges 



For four dimensional black hole with four charges we have c—0, n — 2,b — 1, m — A and 
1 = 1. f(u) is given by 

/ = (1 - u)[l + (1 + fci + k 2 + k 3 + k 4 )u+ 
(1 + k\ + k 2 + k 3 + k 4 + kik 2 + hk 3 + k\k 4 + k 2 k 3 + k 2 k 4 + k 3 k 4 )u 2 - kik 2 k 3 k 4 u z ], 

(D.l) 

and so we have 

F(u) = l + (l + k 1 + k 2 + k 3 + k 4 )u+ 

(1 + ki + k 2 + k 3 + k 4 + k\k 2 + k±k 3 + k x k 4 + k 2 k 3 + k 2 k 4 + k 3 k 4 )u 2 — k 4 k 2 k 3 k 4 u 3 . 

(D.2) 

which is a cubic polynomial. 
From (A. 3) we can write 

Pi(u) = b i0 + c i0 u, (D.3) 
where are constant parameters, cio is given by 

2*4 1 + k ± r k 2 k 3 k 4 

cio = — 610 ^— [T—r b 20 ~ — tt&so - — TT&40 (D.4) 

3 3 1 + k 2 1 + k 3 1 + k 4 

and other c«o can be obtained by permuting (1234) in the above expression. From (A. 5) 
we obtain 

si(u) — hiu + (h 2 + h 3 u) log(l - u/ui) + (h 4 + h 5 u) log(l - u/u 2 ) 

(D.5) 

+ (h 6 + h 7 u) log(l - u/u 3 ) 

Since F(u) is a polynomial of degree 3 the zeroes of F(u), u iy i = 1, 2, 3 have complicated 
expressions. Because of the complicated, which in turn makes explicit expressions of 
parameters hi to h 7 too long to include here. Once again proceeding in a similar way we 
obtain the R-charge conductivities in terms of k^s. 



E Seven dimensional black hole with two charges 



In the case of seven dimensional black hole with two charges we have c = — 1, n = 3, 
6 = 4 and m — 3. 

/ = (1 - u)[l + u+ (1 + h + k 2 + hk^u 2 - fcxAfeu 3 ]. (E.l) 
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So F{u) is given by 

F(u) = 1 + u + (1 + ki + k 2 + A^m 2 - fcifeu 3 , 



(E.2) 



which is, once again, cubic. 

Since / = 2 from (A. 3) we get 

Pi(u) = b i0 + c i0 u 2 . (E.3) 

From (A. 5) we get 

Sl ( u ) = hiu + h 2 u 2 + (h 3 + h A u 2 ) log(l - u/ui) + (g 5 + 5f 6 M 2 )[log(l - u/u 2 ) 

1-^.4) 

+ (97 + 98U )log(l -u/u 3 ). 

F(u) being of order 3, expressions of zeroes i^'s are more involved which makes the 
expressions of hi to h$ too long and complicated to write down explicitly. We have 
analyzed this case and obtain similar conclusion that as critical surfaces are approached 
diffusion coefficients vanish. 
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